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It remains an open question to what extent many of the astronomical sources of intense bursts 
of electromagnetic radiation are also strong emitters of non-photon messengers, in particular grav- 
itational waves (GWs) and high-energy neutrinos (HENs). Such emission would provide unique 
insights into the physics of the bursts; moreover some suspected classes, e.g. choked gamma-ray 
bursts, may in fact only be identifiable via these alternative channels. Here we explore the reach 
of current and planned experiments to address this question. We derive constraints on the rate of 
GW and HEN bursts per Milky Way equivalent (MWE) galaxy based on independent observations 
by the initial LIGO and Virgo GW detectors and the partially completed IceCube (40-string) HEN 
detector. We take into account the blue-luminosity-weighted distribution of nearby galaxies, assum- 
ing that source distribution follows the blue-luminosity distribution. We then estimate the reach of 
joint GW-I-HEN searches using advanced GW detectors and the completed km'' IceCube detector to 
probe the joint parameter space. We show that searches undertaken by advanced detectors will be 
capable of detecting, constraining or excluding, several existing models with one year of observation. 



Gravitational-vifave (GW) astronomy, as well as high- 
energy neutrino (HEN) observations are entering a new 
and promising era with newly constructed detectors, pro- 
viding unprecedented opportunities to observe these as- 
trophysical messengers, opening new windows onto the 
universe. GW observatories [1-4] are being built and 
upgraded to second generation detectors. Several HEN 
detectors [5, 6] have reached their design sensitivities and 
will be further upgraded in the near future [7]. 

GWs and HENs can originate from a number of 
common sources. Plausible sources include gamma-ray 
bursts (GRBs), core-collapse supernovae (CCSNe), soft 
gamma repeaters as well as microquasars [8]. For a 
joint GW-I-HEN analysis, the most interesting sources 
are those which have little or no detected electromag- 
netic (EM) emission. With the near-completion of the 
first searches for multi- messenger GW-fHEN sources [9], 
it is important to examine the projected science reach 
of such searches, as well as how it relates' to indepen- 
dent GW and HEN measurements. This can support 
and guide the theoretical work necessary to gain a better 
understanding of future observational results. 

In this paper we interpret and combine previously pub- 
lished and independent GW and HEN observational re- 
sults, to derive the first joint constraints on the rates of 
GW-I-HEN sources. We first discuss constraints from in- 
dividual HEN and GW searches, and then combine these 
to derive upper limits on GW-I-HEN sources. We finally 
estimate projected constraints on GW-I-HEN sources 
with future detectors and joint GW-I-HEN searches. 

Upper limits from neutrino observations — Abbasi et 
al. [10] searched for transient point sources with the par- 
tially constructed IceCube detector in its 40-string con- 
figuration (hereafter IceCube-40) for over 1 year. The 
search covered the northern sky and emission time-scales 











•-■ 
























— 










■•■ ^— 












— # =; 

— # =; 
— #>= 


3 





= 10 ^ \ ^ ^3 

c 10 10 10 10 

n [#at10Mpc] 



FIG. 1. (Color online) Top: fraction of neutrino-emitting 
sources within 1 Gpc which would be detected with 1, 2, 3, 
or > 3 neutrinos, as a function of Mhen (the mean number of 
detected neutrinos from a source at 10 Mpc) for a detector 
with northern sky coverage (e.g. IceCube). Only sources are 
considered that emit neutrinos towards the Earth. Bottom: 
source population upper limit -R"^ as a function of riHEN, as- 
suming a beaming factor of = 14, and considering only the 
northern sky. 



ranging from seconds to months; no evidence for transient 
sources was found. The number of neutrinos needed for 
detection on each time-scale was determined. Using here 
the conservatively large time window of 500 s for HEN 
emission from GRBs by Baret et al. [11], three spatially 
coincident neutrinos within this time window would have 
been sufficient for a (5cr) discovery. (We note that even 
with the higher event rate of IceCube-86, three coincident 
neutrinos will remain a highly unlikely outcome from the 
background). We therefore estimate the source popula- 
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tion upper limit as the maximum source rate that has 
< 90% probability to result in at least one occurrence 
of > 3 coincident neutrinos in a time window of 500 s 
during a 1-year measurement. 

We model the source population as following the blue- 
luminosity distribution of galaxies [12]: (i) for up to 40 
Mpc, we take the blue-luminosity distribution given in 
the Gravitational Wave Galaxy Catalog [13] (we note 
that any incompleteness in the galaxy catalog makes our 
upper-limits conservative.); (ii) for larger distances (up to 
1 Gpc) we adopt the homogenous blue-luminosity density 
determined by Blanton et al. [14]; (iii) we assume that 
IceCube is uniformly sensitive to sources in the northern 
sky only, which is a reasonable approximation of the de- 
tector's directional sensitivity [10]. Our upper-limits are 
calculated as a function of tihen) defined as the average 
number of detected HENs from a source at 10 Mpc (e.g. 
[15]). 

The results provided here assume that each source has 
the same intrinsic neutrino brightness (limits based on 
a fixed average brightness are conservative compared to 
those using any other brightness distribution), and ac- 
count for beaming of the HEN emission. For a source 
with intrinsic brightness u^en at distance r, the probabil- 
ity that > 3 neutrinos will be detected from it is 

Pin > 3[r, tIhen) = l-F (2[(10 Mpc/r)^^^™) , (1) 

where r is the source distance, F is the Poisson cumu- 
lative distribution function, and n is the number of de- 
tected neutrinos from the source. Therefore for galaxy 
i with blue luminosity at a distance r^, the average 
number Ni of sources which are discovered (i.e. have > 3 
detected neutrinos) will be 

iV,(i?,T) = p{n > 3[r„nHEN) • R/fb ■ T ■ i^'V^B™, (2) 

where R is the source rate [number of sources per year 
per Milky Way equivalent (MWE) galaxy (w.r.t. blue 
luminosity)], fb is the HEN beaming factor of the source, 
T is the duration of the measurement (~ 1 year [10]), and 
L^W ig ^i^e bi^jg luminosity of the Milky Way. The 90% 
confidence source population upper limit R^^ will be the 
upper limit that satisfies 2.3 > iV:,(i?"^ T), i.e. 

n Of T MW 

i?"^(nHBN) = — -J-y (3) 

where the sum is over all galaxies with declination 5i > 0. 
For r > 40 Mpc where we consider a homogeneous matter 
distribution, the summation is substituted with an inte- 
gral. Figure 1 (top) shows the fraction of HEN sources as 
a function of tihen- In the lower plot, population upper 
limits for HEN sources are shown, taking into account 
the sources' HEN beaming factor /{,. As mildly relativis- 
tic jets from CCSNe and low-luminosity (LL) GRBs are 
expected to make up a significant portion of HEN sources 
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FIG. 2. (Color online) Source population upper limits as 
functions of the sources' GW emission in isotropic-equivalent 
energy E^."" . Dashed red: observational limits with ini- 
tial LIGO-GEO- Virgo [18]. Solid blue: projected limits for 
the advanced LIGO-GEO-Virgo GW detectors in the event of 
non-detection. 

of interest [8, 15, 16], we adopt fb = 14 corresponding to 
LL-GRB beaming factor obtained in [17]. 

Upper limits from gravitational-wave observations — 
We use the limits obtained by the latest GW all-sky burst 
search by Abadie et al. [18]. We consider their result for 
sine-Gaussian GW waveform in the sensitive band of the 
GW detectors (- 150 Hz). 

Abadie et al. report no detection using the 
initial LIGO-GEO-Virgo detectors [1-3], and set 
a frequentist 90% confidence upper limit of w 
0.5 X (lO-2M0cV£;^t^)3/7year/Mpc3, or 2.0 detectable 
events per year, on the population of the considered GW 
bursts. Here we interpret this result through introducing 
a GW horizon distance D''-^' {E'^;^) , within which any GW 
bursts with energy would have been greater than the 
loudest background event of the measurement. This is a 
reasonable approximation of the detection efficiency for 
the obtained efficiencies in [18]. The above result yields 
L»™(£;*t?) = 7.8- (i;^tv/10"^^-^o c^)^^^ Mpc. We thus 
derive a galaxy-based GW source population upper limit 
as a function of El^°, using the blue-luminosity- weighted 
distribution of galaxies as described in (i)-(ii) above: 

onrMW i 

V r ^ ' year 

Here, we assumed that each GW source emits the same 
amount of GW energy (limits based on uniform GW 
emission are conservative compared to those using any 
other emission distribution). We estimate the achievable 
population upper limit for the advanced LIGO- Virgo GW 
detector network by assuming a ~ lOx increase in sensi- 
tivity compared to the network of initial detectors, with 
similar measurement duration. Results are shown in Fig- 
ure 2. 
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FIG. 3. (Color online) GW+HEN source population upper 
limits based on the statistical combination of independent 
GW and HEN measurements. Top: observational results 
for measurements with the initial LIGO-GEO- Virgo GW de- 
tectors [18] and the IceCube-40 HEN detector [10]. Bot- 
tom: projected results for 1-year observations with advanced 
LIGO-Virgo and IceCube-86. The limits shown assume a 
HEN beaming factor of 14. Horizontal lines: expected HEN 
rate from the Waxman-Bahcall [19] (solid) and Ando-Beacom 
[16] (dashed line) models, scaled to the IceCube-40 (top) and 
IceCube-86 (bottom) detector configurations. 



Joint GW+HEN population upper limits — Individual 
GW and HEN observations can be combined to deter- 
mine a GW-I-HEN source population upper limit in the 
Eq^^-u^ien parameter space. In Figure 3 (top) we provide 
GW-I-HEN population upper limits based on the statis- 
tical combination of current observational results from 
independent GW and HEN measurements. We obtain a 
joint observational upper limit by considering that, on av- 
erage, less than 2.3 GW-I-HEN bursts occur within 
or have > 3 detected HENs per year (this is a > 90% 
confidence upper limit, since the GW and HEN mea- 
surements were longer than one year). The observational 
GW-I-HEN upper limit for a source population propor- 
tional to the blue-luminosity-weighted galaxy distribu- 



tion will therefore be 
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X]{ri>_D"w,5i>0} — Sj?"!, ".HEN)-t'B + S{ri<_Df-w} ^ 

(5) 

Figure 3 (top) also compares the obtained upper limits to 
two HEN emission models ([16, 19]). As the theoretical 
estimates in [16, 19] are provided for km'' scale detectors, 
we convert them to estimates for IceCube-40 by conser- 
vatively assuming that the sensitivity of IceCube-40 is 
one-third that of IceCube-86. 

Similarly to the above observational results, we also 
calculate the projected GW+HEN population upper lim- 
its based on the statistical combination of projected re- 
sults from independent, 1 year long measurements with 
advanced LIGO-GEO- Virgo and IceCube-86. Results are 
shown in Figure 3 (bottom). 

We now estimate the projected population upper limits 
for GW+HEN sources obtainable with a jomt GW+HEN 
search, considering a 1-year measurement with the ad- 
vanced LIGO-Virgo and IceCube-86 detectors. We con- 
sider an event candidate to be the coincidence of 1 
GW trigger and 1 HEN. While we might detect more 
than 1 HEN from some sources, the fraction of such 
sources is small (see Figure 1), therefore we conserva- 
tively omit multi-HEN sources. For the joint search we 
define a horizon distance D°*"™(£^^^), such that a joint 
GW+HEN event with 1 detected HEN and GW energy 
E^f^, within E)™^™ would be more significant than the 
(anticipated) loudest background event. We estimate 
£)GWHEN ^Yie same as the exclusion distance of the 

externally triggered search for GW bursts by Abbott et 
al. [20], who obtained a median exclusion distance of 
L> ~ 12 Mpc • (£;^t^/10-2Moc2)i/2 with GW emission in 
the most sensitive band of LIGO-Virgo. Such comparison 
to externally triggered GW searches is a reasonable ap- 
proximation if the joint search has 0(1) chance overlaps 
of background GW and HEN events (which can be con- 
trolled by adjusting the event selection threshold). For 
the joint GW+HEN search the estimated source popula- 
tion upper limit R^'^ will be 
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(6) 



The estimated population upper limits for a GW+HEN 
search are shown in Figure 4 for advanced LIGO-Virgo 
and IceCube-86. 

Discussion. — To compare our results to emission mod- 
els, we consider SNe with mildly relativistic jets (SNe 
with jets) as promising GW+HEN emitters, whose ob- 
servational rate upper limit has been about < 1% of the 
total SN Ib/c rate [21] (- 3 x IQ-^/year/MWE galaxy). 
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joint search is > 2.4 times greater 
individual GW searches. 
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FIG. 4. (Color online) Projected GW+HEN source popula- 
tion upper limits for a joint analysis of one-year of observa- 
tions with advanced LIGO- Virgo and IceCube-86. Results 
are given as functions of source emission parameters Eq^ 
(GW emission in isotropic-equivalent energy) and tihen (aver- 
age number of detected neutrinos from a source at 10 Mpc). 
Horizontal lines: expected HEN rate from the Waxman- 
Bahcall [19] (solid) and Ando-Beacom [16] (dashed line) mod- 
els. The limits shown assume a HEN beaming factor of 14. 



All-sky population upper limits with IceCube-86 are 
projected to exclude sources at the rates of SNe with 
jets, for riHEN ^ 300 (Fig. 1). This is comparable to the 
Waxman-Bahcall flux [19], which is, in the specific case of 
GRBs, already constrained by HEN searches [22]. Upon 
non-detection, advanced GW detectors are projected to 
exclude sources at current limits on rates of SN with jets, 
for > 5 X IO-^MqC^ (Fig. 3, bottom panel), exclud- 
ing the suspended accretion model [23] , and constraining 
the accretion disk fragmentation model [24]. 

The above estimates assume that the fraction of SNe 
with jets is probed by radio observations. It has been 
proposed, however, that mildly relativistic jets may be 
much more common, but completely choked (bright in 
neutrinos, dark in gamma-rays and radio) [15, 16]. The 
nearby core-collapse supernova rate is high enough to 
allow testing these models soon. Upon non-detection, 
IceCube-86 is projected (Fig. 1) to exclude sources at SN 
rates with tIhen ^ 12, a level comparable to the emission 
expected from SN jets by Ando and Beacom (tihen ~ 
10; [16]), or emission through reverse shocks in mildly 
relativistic jets (tihen ^ 7; [15]). Moreover, as evident 
from Fig. 2, advanced GW detectors are projected to 
exclude sources at SN rates with average GW emission 
of ^ 2 X 10~'^Mqc'^ , excluding suspended accretion 
as well as accretion disk fragmentation. 

We obtain projected population constraints with a 
joint GW-I-HEN search (Fig. 4) that can be more re- 
strictive in some regions of the parameter space than in- 
dividual searches if the GW horizon distance D™ of the 
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